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ABSTRACT 



This work presents high spectral resolution observations of the [C n] line at 158 /im, one of the major cooling lines of the interstellar 
medium, taken with the HIFI heterodyne spectrometer on the Herschel satellite. In BCLMP 691, an H n region far north (3.3 kpc) in the 
disk of M 33, the [C n] and CO line profiles show similar velocities within 0.5 km s -1 , while the H i line velocities are systematically 
shifted towards lower rotation velocities by ~ 5 km s -1 . Observed at the same 12" angular resolution, the [C n] lines are broader than 
those of CO by about 50% but narrower than the H i lines. The [C n] line intensities also follow those of CO much better than those 
^5, of Hi. A weak shoulder on the [Cn] line suggests a marginal detection of the [ 13 Cn] line, insufficient to constrain the [Cn] optical 

depth. The velocity coincidence of the CO and [C n] lines and the morphology at optical/UV wavelengths indicate that the emission 
is coming from a molecular cloud behind the H n region. The relative strength of [C n] with respect to the FIR continuum emission 
is comparable to that observed in the Magellanic Clouds on similar linear scales but the CO emission relative to [C n] is stronger 
in M 33. The [Cn] line to far-infrared continuum ratio suggests a photoelectric heating efficiency of 1.1%. The data, together with 
published models indicate a UV field Go ~ 100 in units of the solar neighborhood value, a gas density n H ~ 1000 cirT 3 , and a gas 



temperature T ~ 200 K. Adopting these values, we estimate the C + column density to be N c + ~ 1.3 x 10 17 cirT 2 . The [Cn] emission 
comes predominantly from the warm neutral region between the H n region and the cool molecular cloud behind it. From published 
abundances, the inferred C + column corresponds to a hydrogen column density of N H ~ 2 x 10 21 cm" 2 . The CO observations suggest 
that Nh = 2Nh 2 ~ 3.2 x 10 21 cirT 2 and 21cm measurements, also at 12" resolution, yield Nm ~ 1.2 x 10 21 cirT 2 within the [Cn] 
■ velocity range. Thus, some H 2 not detected in CO must be present, in agreement with earlier findings based on the SPIRE 250 - 500 

f^*) 1 yum emission. 

Key words. Galaxies: Individual: M 33 - Galaxies: Local Group - Galaxies: evolution - Galaxies: ISM - ISM: Clouds - Stars: 
Formation 
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1. Introduction [Cn] emission can originate in otherwise neutral gas. A se- 

rrt 1 ries of [C n] observations in the Local Group galaxy M 33 

The [Cn] line at 157.741 /an is the strongest emission line with the fflFI (Heterodyne Instrument for Far Infrared astron- 

from galaxies, representing close to one per cent of the bolo- omy) instrument on th e Hers chel Space Obse rvatory (iPilbrattl 

metric luminosity (ICrawford et al.| [L985J; I Stacey etal. | [19911 Ide Graauw et all l20lH: iRoelfsema et al1l20llh is the es- 

Malhotraetal.||2001|) and as such is detectable in very dis- sential comDonent of me Qnen Time Kev Project HerM33ES 



iviainoira ei ai.| |zuyi r aim as sucn is uereciaoie in very uis- sen tial component o f the Open Time Key Project HerM33ES 
tant objects <|LoebJ [19931 | Maiolino et al. | [2005|). However, the ([Kramer et al.ll20lol) . A first article dMookeriea et al.ll20ll pre- 



question of the origin of the [Cn] emission from galaxies sented a pACS map of me [Cn] and [Ql] emission and a H IFI 

dHeiles| [1994J) remains open to this day. Possible sources are [c n] spectrum of the H n re g ion BCLMP 302. This work is the 

the diffuse warm ionized medium, H n regions, and mostly neu- first of a series preS enting extended high spectral resolution ob- 

tral photo-dissociation regions (PDRs). Because carbon is ion- servations of the [Cn] line in M 33 
ized more easily than hydrogen (11.26eV rather than 13.6eV), 

Despite its importance, only a few velocity -resolved obser- 

* Herschel is an ESA space observatory with science instruments vations have been made of the [C n] line (e.g. |Boreiko & Betj 

provided by European-led Principal Investigator consortia and with im- 119911) . The goal of such observations is not only to obtain line in- 

portant participation from NASA. tensities but also velocity emission profiles to enable a compari- 
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son of the velocities and velocity widths of the [C n] -emitting gas 
with those of the lines due to CO and HI (and Ha in the future). 
The [C n] line emission provides information on the heating and 
cooling processes of the gas and is a rather unique probe of the 
photon-dominated region (PDR) between the cool CO-emitting 
molecular gas and the gas ionized by massive stars. This region 
is likely to be more massive and extended in sub solar-metallicity 
galaxies and may represent a significant but hitherto unmeasured 
gas mass. Hy per-fine component s of the [ 13 Cn] line bracket the 
main isotope (ICooksv et all 19861) and with HIFI we are, in prin- 
ciple, able to separate the lines allowing a direct determination 
of the optical depth. 
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Fig. 1. North-east ern part of M33 in the IRAC 8 ^m band 
dVerlev et al.ll2007l) . The upper (small) rectangle shows the re- 
gion presented in Fig. [2] and the small triangles forming an 
'X' in and around the rectangle indicate the individual posi- 
tions of the [Cn] HIFI observations. The large lower rectangle 
shows the region aroun d BCLMP302 observed with PACS by 
iMookeriea et al.l d201 ll) . The center of M 33 is at the position of 
the triangle to the lower right of the image of M 33. The giant 
H ii region NGC 604 can be seen as the very bright region to the 
left around Dec 30:47 and is marked with a x. 



Recently, lLanger et al.l (1201 Ol) found a population of diffuse 
interstellar clouds (Ay < 1.3 mag) in the Milky Way which 
show [Cn] and Hi emission, but no detectable CO emission. 
The observed [C n] emission is stronger than expected for dif- 
fuse atomic clouds. The idea of CO-dark H2 is not ne w and (a) 
some is definitely found ar ound CO-bright clouds (e.g. lPak et al.l 
119981: IWolfire et aDl2010h but pos sibly also (b) mixed with the 
H 1 (e.g. iPapadopoulos et al.l 20021) or (c) i n the outermost parts 
of galactic disks (e.g. Pfenniger et al.ll994l) . lLanger et al.l (l2010h 



attribute the [Cn] excess emission to diffuse and warm (T^ 100 
K) molecular (H2) clouds. While the linear resolution of these 
observations is excellent, the Galactic plane is optically thick at 
short wavelengths such as the Ha line commonly used to trace 
Hn regions. Thus, some of the [Cn] emission they detect may 
in fact be due to ionized rather than neutral gas. The observa- 
tions presented here allow us to address the same questions with 
~ 50pc resolution but with the advantage of much clearer views 
of the region and of the sources of emission at all wavelengths. 

Previous observations have established that [Cn] emission 
is a good tracer of star formation regions. Among Local Group 
galaxies, there is an excellent correspondence between \C ill 
and Far-IR emission in b oth Magellanic Clouds dPoglitsch et al.l 
119951: llsrael et al.lll996t llsrael & Malonevll201 lb and between 
[C n] and both Ha and mid-IR (24 /jm ) emissio n in the spiral 
arms of M 31 (Ro driguez -Fernandez et al.1l2006l) . The fact that 
the [C 11] line represents a similar fraction of the global FIR lumi- 
nosity of spiral galaxies shows that this also holds in a statistical 
way on galactic scales. 

As illustrated in Fig. 1, the HerM33es project observes [C 11] 
in a long strip along the major axis. In this work, we investigate 
the relationship between [C 11] FIR, CO, and H 1 emission at the 
scale of individual clou ds (~50 pc) in an area centered on the 
H 11 region BCLMP 6 91 dBoulesteix et al.ll 19741) and compare it 
with the results from IMookeriea et al.l (1201 ll) on BCLMP 302, 
both part of the strip. In particular, we estimate the amount of 
gas traced by the [C n] emission and and relate that to the CO 
and HI emission. 

2. Observations and data reduction 

The HIFI Wide Band Spectrometer horizontal and vertical 
polarization data were processed to level 2 using the stan- 
dard pipeline of HIPE version 4 and then exported to CLASS 
(see http://www.iram.fr/IRAMFR/GILDAS). The obsids for the 
BCLMP691 observations are 1342213155 and 1342213156 
which are load-chop On-The-Fly maps consisting of 8280 in- 
dividual spectra each. The maps are perpendicular strips with 
position angles of 22.5 and 1 12.5 degrees with close spacing be- 
tween individual spectra. All spectra were smoothed to 8 MHz 
(1.26 km s~') resolution and a zero-order baseline (i.e. a con- 
tinuum level) was subtracted excluding velocities from -280 - 
-255 km s , corresponding to the expected line velocity. The 
baseline was calculated over a region extending from -345 to 
-170 km s _1 . All velocities are in the LSR frame and the sys- 
temic velocity of M 33 is about -170 km s _1 . 

A number of spectra (roughly 1/3) suffer from the HIFI 
HEB (hot electron bolometer) standing wave pattern at about 
290 MHz. Various Fourier filtering schemes were applied but it 
proved very difficult to reliably eliminate the standing waves. 
In the end, the 11298 good spectra were kept and the others 
dropped. 

At 5" intervals along the cross, the spectra were convolved 
with a 5" gaussian in order to obtain a single high quality spec- 
trum at each position. The resulting beam size is 12.3x1 1 .4" be- 
cause the convolution only stretches the beam in the direction 
where there are spectra (i.e. along the lines of the cross). In the 
center, the beam is more circular because both arms of the cross 
are included in the convolution. Four "extra" spectra were gen- 
erated near the center at (±5", ±5") with respect to the cross (see 
Fig. |2). Although no pointings correspond exactly to these po- 
sitions, this is possible because the neighboring beams overlap 
considerably in this region so these points are sampled on both 
arms of the cross. 
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Velocities, line widths, and line intensities were mea- 
sured from gaussian fits to the 17 detected positions on 
the cross plus the 4 "extra" positions. Table 1 gives the 
positions of the spectra with the fit results for the de- 
tected positions and the rms noise level for the others. 
The uncertainties come from the gaussian fits using the 
CLASS fitting routine. As per the current HIFI users man- 
ual http : //herschel . esac . esa . int/Docs/HIFI/html/ 
chQ5sQ5 .html#beam-in£ormation, we adopt a main beam 
efficiency of 0.693 and present all data on the main beam 
temperature scale. 

3. Analysis of BCLMP 691 [Cn] line intensities 

3.1. The H ii region complex BCLMP 691 

Both BCLMP 691 (IK 60, Aller 3, Searle 15) and BCLMP 302 
(als known as IK 53, Aller 43, Searle 11) are among the 
brighter nebulae in the northern half of M 3 3 with 1.4 GHz 
flux d e nsities of 4.6 and 6. 1 mJy res pectively ( dBoulesteix et al.l 
119741) . dlsrael & van der Kruiill 19741) ). Figure 1 shows the north- 
ern part of M 33 with the region observed around BCLMP 691 
(and BCLMP 302) indicated. Zooming in on BCLMP 691, 
Figure [2] shows the fCnj det ections superposed on a 24 yon 
image (TTabatabaei et al.ll2007l) of the region. BCLMP 691 and 
BCLMP 302 are similar in many aspects. Both are about 75 
pc in diameter with large-scale r.m.s electron densities of < 
>'/ 2 of 7.4 and 6.3 cirT 3 and excitation parameters of 



Rstr 



omgren rl e 



180 and 200 pc cm , corresponding 



to excitation by the equival ent of 3 and 4 05 exciting stars 
dlsrael & van der Kruitlll974l) . More significant differences be- 
tween BCLM P 691, the HII region studied in this paper, and 
BCLMP 302 ( dMookeriea et al.ll201 lT) t are tied to their location. 
BCLMP 691 is a major Hn region 3.3 kpc from the center, and 
relatively isolated, whereas BCLMP 302 is a major star-forming 
region embedded in the northern spiral arm at 2.1 kpc from the 
center of M 33. The Hi and CO(/=2-l) emission towards the 
latter is stronger than towards BCLMP 6 91. BCLMP 691 ha s 
a metallicity 12+log(0/H) = 8.42 ± 0.06 dMagrini et alJl2010h . 
i.e. an oxygen abundance relative to H of 2.6 x 10~ 4 , about 
half the metallicity of the solar neighborhood. This is some- 
what below t he average of the m ore central Hn regions such 
as NGC 595 dMagrini et al.ll2010h . BCLMP 302 is irradiated by 
the combined young and old stellar population of the northern 
spiral arm, but in the vicinity of BCLMP 69 1 , the radiation field 
from the old stellar population is weak. Nevertheless, the differ- 
ence in ionization between the two regions is remarkable: the 
isolated complex BCLMP 691 has a much higher degree of ion- 
ization than BCLMP 302, with ratios Ne ++ /Ne + = 1.25 (0.19) 
and S 3+ /S ++ = 0.15 (0. 04) where the valu es in parentheses are 
those of BCLMP 302 dRubin et al]l2008l) . The ionization dif- 
ference may be attributable to the lower metallicities at greater 
radii causing the exciting stars to have harder radiatio n fields. 
As can be seen in the appendix of iGratier et al.l (1201 ll) . the ge- 
ometry of the H ii region - molecular cloud interfa ce differs be- 
tween BCLMP 302 and 691 (clouds 256 and 299 of IGratier et al] 
1201 ll) . It is clearly of interest to identify the corresponding ef- 
fects on the [C n] emission. 

3.2. The [C \\] to FIR and the [C \\] to CO ratios 

The large-scale [Cn] to FIR continuum flux ratio in nor- 
mal galaxies i s usua ll y below 1%, typi c ally 0.3 - 0.5% 
dCrawford et al.l 119851: fstacev et al.l 1 199 it iBraine & Hughesl 



119991: lMalhotraetaLll200ll) . Here, FIR refers to the dust emis- 
sion in the wavelength range 42-122/mi (IRAS definition). This 
is about half of al l dust emission (i.e. FIR » 1/2 TIR, see 
iDale & Heloull2002l) . At the scale of a spiral arm in M 31, using 
ISO L WS data at 300pc resolution, iRodrfguez-Fernandez et al.l 
J2006I) find a considerably higher ratio (I[c 111/FIR42 -122 ~ 2%). 
The cl oses t to M 33 in metallicity is the LMC, where llsrael et al.l 
d 19961) and llsrael & Malonevl d2011l) find I C n/FIR flux ratios be- 
tween of 0.7 and 5.0 %. The [C n] flux can be calculated as 
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Tdv 



(1) 



where J Tdv is the integrated intensity in the standard units of 
K km s (the 10 s converts from K km s to K cm s ). 

For BCLMP 691 we find a ratio of I [C ii]/FIR 4 2-i22 ~ 1.1% 
and a ratio of [Cn] flux to total dust emission of 0.54%. 
The Fciii/IR ratio is i ndicative of the gas heating efficiency 
dBakes & Tielenslll994l) and the value for BCLMP 691 i s very 
close to the average ratio obtained by iMookeriea et alJ d201 ll) 
for BCLMP 302 (their region A). In this discussion we focus on 
the center of the H 11 region, our (0,0) position. The FIR and to- 
tal dust luminosities for this region were calculated by directly 
summing the emission in these wavebands using the 5.8, 24, and 
70 yt/m emission from Spitzer and the PACS and SPIRE data ob- 
tained by the HerM33es project. PAH emission has been left out 
of the total dust flux. In both M 33 fields observed, these ratios 
are thus very simila r to each other and to th ose obtained in the 
LMC and the SMC dlsrael & Malonevll201 ll) . 

We may also compare the [Cn] to CO flux ratios for the 
various environments. Unfortunately, we do not have 7=1-0 
12 CO measurements at 12" resolution so we have convolved our 
CO(2-l) observations to the 25" resolution of the CO(1-0). At 
this lOOpc resolution, the CO(2-1 )/(1-0) line ratio is .72, very 
close to the average value given in IGratier et al.l d2010l) . 

For BCLMP 691, a comparison of [C n] and CO lines can be 
found in Tab. 1 and the spectra are superposed in Fig. [3] As for 
[Cn], the CO flux can be calculated as I co = 10 5 J Tdv = 
1.57 x 10~ 9 J Tco(i-Q)dv where J Tdv is in Kkms -1 as usual. 
Assuming the CO(2-1)/(1-0) line ratio is also 0.72 at 12" res- 
olution, we find a I[Cn]/Ico(i 0) nux ratio of about 8000 for 

BCLMP 691, slightly lower than the BCLMP 302 peak value 
(Mookerjea et al). Although this is in excess of the ratios of 
about 1500 found in normal spiral galaxies, and even of ratios 
of 6000 such as found i n starburst galaxies dStacev et al.ll 199 it 
IBraine & Hughesl fl999). it is considerably lower than the cor- 
responding ratios found in the LMC, which (with one excep- 
tion, N 159W) range from 13000 to more than 100000 with a 
mean of 60000 for the individual star-formin g clouds observed , 
and a ratio of 23000 for the LMC as a whole dlsrael et al.lll996t 
IIsrael&Malonevll2011l) . 

In general, because of its (indirect) dependence on UV pho- 
tons, [Cn] emission is linked to star formation. GMCs with- 
out H 11 regions have much weaker [C 11] emission for compa- 
rable levels of CO em ission, such as N 159S in the LMC with a 
I[Cii]/Ico ratio of 400 dlsrael et al.lll996l) . In comparison with the 
LMC clouds, the two objects in M 33 have rather low I[Cn]/Ico 
ratios although I[Cnj/FIR ratios are practically identical. Since 
the latter are a measure of the heating efficiency of the clouds, 
the conclusion must be that otherwise comparable star-forming 
complexes are much richer in CO in M 33 than in the LMC. 
Since there is no difference in this respect between the two M 33 
objects, the degree of ionization, i.e. the hardness of the radiation 
field, is not a factor of significance. 
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3.3. Carbon column density 
3.3.1. A lower limit to N c + 

The [C n] line is an important probe of the physical conditions 
of the ISM. In particular, it provides information on regions not 
probed by the CO line and possibly not probed by the H i line. 
The [C n] line is thermalized at a density of about 3000 cirT 3 in 
neutral gas and has an upper level energy of 91 K. Thus, a mini- 
mum [C n] column density can be calculated assuming that den- 
sities » 3000 cm -3 and temperatures » 9 IK. Lower densities 
and/or lower temperatures would decrease the [C n] emission per 
H-atom. A lower limit to the hydrogen colu mn density can be 
dete rmined by assuming a car bon abundance. Ilsrael et al.ld 19961) 
and Ilsrael & Malonevl (1201 ll) performed comparable studies at 
very similar linear resolutions in the brightest H n region com- 
plexes of the Large Magellanic Cloud, which has a metallicity 
similar to that of M 33. These studies therefore provide an excel- 
lent basis for comparison, even though they did not measure the 
actual [C n] l ine profiles. Fol l owing the reasoning and equa tions 
presented bv lCrawford et al.1 (1 19851) and Ilsrael et aD (Il996h - 



^ > 6.34x 10 20 / |C „]cm- 2 



(2) 



with /[en] defined in Eq. (1). 

The [C n] line intensity at the peak of 9.9 K km s _1 is equiv- 
alent to 7.0 x 10~ 5 erg cm _2 s _1 sr _1 , or a column density of 
Nc* > 4.4 x 10 16 cirT 2 . We emphasize once more that this is 
a lower limit because of the high-density, high-temperature as- 
sumption underlying these numbers. We have also assumed that 
all of the ionized carbon was singly ionized - this latter hypoth- 
esis is reasonable for the neutral medium (i.e. where H is not 
ionized). 

That carbon is only singly ionized is not a good approxi- 
mation for regions where Hn is the dominant phase. The ion- 
ization potential for doubly-ionized carbon is 24.38 eV, consid- 
erably below that of doubly-ionized neon (40.96 eV) or triply 
ionized sulphur (34.83 eV). Given the high degree of ionization 
characterizing BCLMP 69 1 (Section 3 . 1 ), in the H n region itself 
most carbon will be in the higher ionization states. Therefore, in 
the following, we will compare the derived carbon column den- 
sity with that of the neutral gas. By the same token, we expect 
little of the [C n] emission to be contributed by the ionized gas 
in the H n region, which would have caused us to overestimate 
the amount of emission from the neutral gas. Rather, we may 
assume that essenti ally all of the fC i ll is in deed associated with 
the Hi and FL;. As iMookerjea et all d201 lb have estimated that 
emission from the H n region BCLMP 302 might account for 20 
- 30% of the [C n] emission, this is another potential distinction 
between that complex and the one studied here. 
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Fig. 3. [Cn], CO, and Hi spectra of all detected positions at 
12" resolution. Offsets in arcsec are indicated in the upper right 
corner of each panel. The order follows that of Table 1 from 
bottom left to upper right. [C n] is in black, CO in green, and H i 
is in blue. Y-scale is main beam antenna temperature in Kelvins 
for [C n] and CO; H i temperatures need to be multiplied by 150 
to obtain true antenna temperatures. 



3.3.2. A more realistic estimate of N c + 

We can get an idea of how much of an underestimate Eq. (2) is 
by looking at Fig. [4] which shows ionized carbon column den- 
sity variations as a function of density and temperature for low 
optical depths. Under these conditions the C + column density is 
proportional to the [C n] line intensity, so that we can use this 
diagnostic diagram to calculate C + column densities when infor- 
mation about the temperature and density of the gas is available. 
Throughout, the ratio between Nc+ and I[Cn] is at least twice that 
of Eq. (2), except in the very high temperature, very high density 
regime (the violet region at the upper right). 

It is unlikely that the gas (or the dust) in BCLMP 691 
has a temperature T» 9 IK and density n » 3000 cm -3 on 



the (linear) scales we are sensitive to here. As no t ed be fore, 
these scales are almost identical to those Ilsrael et al.l (1 19961) and 
Ilsrael & Malonevl (120111) were sensitive to in the LMC. In the 
LMC regions N 159, N 160, and Nil they estimate Go to vary 
between 25 and 450. Assuming the flux impinging on a dust 
cloud is I star - Anldust, where the An converts the ergcm^s -1 
sr~' to erg cm _2 s _1 , that hal f of the stellar flux I star is du e to 
photons with hv < 6eV (e.g. iTielens & Hollenbachlll985l) . we 
can es timate the FUV (6 < hv < 13.6eV) flux expressed in units 
of the|HabinJ dl968j) field to be G = 27r^„ Jf /1.6 x 10 3 . 

Direct integration of the dust emission from 5 to 610 /urn (the 
edge of the SPIRE 500 fim band) yields a peak flux of 0.013 
ergcm _2 s _1 sr~\ corresponding to Go = 51, well within the 
range found to apply to the LMC clouds. Since the core of the 
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Fig. 4. Plot of the log of the C + column density as a function 
of gas density and temperature for the [C n] intensity at the (0,0) 
position of I[Cn] = 7.0 x 10~ 5 erg cm _2 s _1 sr . The bottom of 
the color scale (black) shows the high-temperature high density 
limit from Eq. 2 which yields log(/Vc+) > 16.65 - in the range 
of parameter space presented here, the limit is not reached. The 
column densities are for optically thin emission. Contours are 
plotted at column densities of log(Nc+) = 1 7, 17 .3, and 17.6. 
We have used Eq. Al from ICrawford et all (1 19851) assuming a 
critical density of 3000 cirT 3 . 



Hn region BCLMP 691 is smaller than the 12" beam (see Ha 
contours in Fig. [2] the central Go is likely to be somewhat higher 
but probably lower than in N 159 or N 160. At these temperatures, 
the dust should emit strongly in the Spitzer 24 and 70yum bands; 
if the grains are in thermal equilibrium, then a temperature of 
about 90 K is suggested by the flux rat io and the gas kinetic 
tempe rature should be somewhat higher dHollenbach & Tielensl 
119991) . The strong 5.8/zm emission shows that some warmer dust 
is also present. In N 159 and N160. llsrael et al.1 (1 19961) estimate 
that the C + column density is 2 - 4 times higher than the limit 
given by Eq. (2). As Go rises, we expect the gas temperature to 
rise, moving C + column densities towards the limit expressed by 
Eq. (2). Hence, for BCLMP 691, our best estimate is a column 
density Nc* ~ 1.3±0.4xl0 17 cm -2 , 3 times the lo wer limit in Eq. 
2). B elow, we confirm this estimate using the iKaufman et alJ 
19991) PDR model. 

3.3.3. Hydrogen column density N w 

There are no direct mea surements of the carb on abundance in 
BCLMP 691. However. lEsteban et all d2009l) estimate carbon 
and oxygen abundances for NGC 595 and NGC 604 and the 
C/O ratios they fi nd are \C]/\Q] = -0.16 and -0.20. They 
present models by ICarigi et all (|2005) that suggest C/O ratios 
of about [C]/[<9] = -0.33. Based on the models and these re- 
cent and high-quality data, we assume [C]/[<9] = -0.25 ± 0.08 
yielding a carbon abundance for BCLMP 691 of 12+log[C]/[H] 
= 8.17 ± 0.1, or a C abundance of 1.48 x 10~ 4 +25% after prop- 
agating the (independent) uncertainties given above for the O 
abundance and the C/O ratio. If we assume that half of the car- 
bon is locked in dust grains, this leads us to a gas-phase C abun- 



dance of 7.4 x 10~ 5 . Therefore, for /V c + ~ 1.3 X 10 17 cirT 2 , the 
H column density is about Nu * 1.8 x 10 21 crrT 2 . If the frac- 
tio n of carbon lock ed in dust is 2/3 (as suggested by Table 2 
of LSofia etal1l2011h , then the gas-phase C abundance would be 
correspondingly lower and the H column density would become 
Nn * 2.6 x 10 21 cnr 2 . 

We can compare this column density with that of the HI 
and the CO. The Hi column density at 12" resolution towards 
the center of BCLMP 691 is N m = 1.7 x 10 21 cirT 2 , of which 
1.2 x 10 21 cirT 2 is within the [Cn] velocity range. Assuming a 
conversion factor 7V(H2)//co = 4x 10 20 cirT 2 /(K km s _1 ), twice 
that of the Galactic disk, the molecular gas column density is 
N„ = 2N(H 2 ) ~ 3.2 x 10 21 crrT 2 . The [Ci] line has not been 
observed in BCLMP 691. At larger scales, the whole Hn re- 
gion emits about 2% of the Ha flux of M 33 and the associated 
molecular cloud 2 99, which is about t he same size as the H n re- 
gion (see table in iGratier et al.ll2.01 lT) . emits about 0.1% of the 
CO emission. 

4. Analysis of BCLMP 691 [C n] line profiles 

4.1. The 13 CII line 

The three hyperfine components of the [ 13 Cn] line have fre- 
quencies equivalent to offsets of -66, +11, and +63kms~', 
respec tively, from the [ 12 Cii] line center at 1 900536.9 + 1.3 
MHz (ICooksv et al.1 11986k IStacev et al.l 1 199 it iBoreiko & Betzl 
[19961) . The ratio of the sum of the three hyperfine [ 13 Cn] tran- 
sitions to the [ 12 Cn] line intensity, coupled with the relative 
abundance ratio, prov ides a measure of the optical depth of the 
[ 12 Cn] emission (e.g. IStacev et al.lfl99ll) . We have searched for 
a possible signature of [ 13 Cn] emission in the inner regions of 
BLMP 691 where the +11 kms -1 component may occur as a 
shoulder present on the [ 12 Cn] line. After fitting a gaussian to 
the sum of the five strongest detections (the 0,0 spectrum and 
the first spectrum on each leg of the cross), we find that the 
gaussian fit is good on the blue side (-280kms _1 ) but leaves 
an excess on the red side, at 10-11 km s _1 from the [Cn] line 
center (see Fig- [5j. This transition, F = 2 — » 1, i s expected to 
be the strongest of the three hyperfine transitions (ICooksv et al.l 
Il986l) . Although the velocity coverage of the spectra includes the 
two other components, they are not detected. The lower panel 
of Fig.|5]shows the residual (spectrum minus gaussian fit) with 
the nominal position of the [ 13 Crx] line. Using 10 km/s windows 
shifted by -66, +11, and +63 km s _1 with respect to the [ 12 Cn] 
line center, we find a nominal summed [ 13 Cn] line intensity of 
0.31+0.23 K km s _1 , i.e. a ratio 97 > /12/13 > 16 with one sigma 
uncertainties. 

4.2. [Cw] line profiles compared to CO and H 1 

As can be seen from both Table 1 and Fig. [6] the CO and [Cn] 
radial velocities determined from gaussian fits agree very well, 
generally to within half of a CO velocity channel. The excep- 
tion is the relatively poor measurement at position (+7.7, +18. 5) 
which is a (3.5cr) detection with a very large uncertainty asso- 
ciated with both the velocity and line width. In BCLMP 302, 
Mookerjea et al. (201 1) find a mean velocity difference between 
CO and [Cn] of 1.3 km s , presumably reflecting physical mo- 
tion (outflow?) of one of the gaseous components giving rise to 
the [Cn] emission. We do not find such a velocity difference 
in BCLMP 691 but that could well be a question of geometry. 
The [C n] profiles presented here are similar in width to the HIFI 
spectrum presented by Mookerjea et al., as are the CO profiles. 
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-300 -250 -200 

Velocity (km/s) 

Fig. 5. [C n] spectrum of the sum of the 5 strongest positions 
of the cross and a gaussian fit (green), showing the shoulder 
at v « -259kms _1 and the velocity where [ I3 Cn] emission 
would be expected. The average CO spectrum of these positions 
is shown as a dashed line. The lower panel shows the [C n] resid- 
ual (spectrum minus gaussian fit). 



The [Cn], Hi, and CO spectra for all positions detected in 
[Cn] are presented in Fig. [3] While Hi is detected everywhere, 
the CO emission peaks in the regions where [C n] was detected. 
The CO and [C n] profiles are similar in shape although the [C n] 
lines are broader. H i, on the other hand, is systematically shifted 
to more positive radial velocities (i.e. lower rotation velocities) 
by about 5 km s and the CO and [Cn] profiles are at the high 
rotation velocity edge of the H i. Only a few of the [C n] spec- 
tra show emission at velocities similar to those of the Hi (e.g. 
low rotation velocity tail of emission at offsets 9.2,-3.8 and 
13.9,-5.7). 

Figure |7]presents a comparison of the position-velocity (PV) 
plots of [Cn], CO, and Hi along the two arms of the cross ob- 
served in [Cn]. In both directions, the CO temperature and CO 
velocity follow those of [C n] closely but the H i profiles seem 
uncorrected in temperature and spread over a much broader 
range in velocities (typically ~ 20 km s~'). Figure [6] shows the 
general agreement in velocities between CO and [C n] in a dif- 
ferent way but it also suggests (lower panel) that the line widths 
mostly fall into two groups. The brighter regions (smaller error 
bars) define a distribution where the CO line width is essentially 
constant while the [C n] line width varies by a factor two. At the 
very least this indicates that the CO-emitting gas is dynamically 
cooler, less affected by the processes that cause large-scale mo- 
tions in the bright [C n] emitting region. The weaker regions, on 
the other hand, show a pattern where both CO and [C n] appear 
to vary linearly with each other. These regions correspond to the 
more diffuse ISM. Looking more closely at the bright regions 
(i.e. ignoring the regions with larger error bars) we see that there 
may be a division into two further groups with [C n] widths of 
~ 10 km s and ~ 15 km s respectively and CO widths of 
~ 7 km s _1 and ~ 8 km s _1 respectively. The broad and strong 
[C n] lines are found along the Hcr-bright extension to the SE. 
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Fig. 6. (top) Comparison of CO and [Cn] velocities, showing 
the good agreement, (bottom) Comparison of CO and [C n] line 
widths, showing that with the exception of some positions with 
large uncertainties, the CO lines are narrower. The velocities, 
line widths, and their uncertainties, have been measured by gaus- 
sian fits to the spectra. 



5. Synthesis and Conclusions 

In order to better locate BCLMP 691 on the diagnostic diagram 
in Fig. HI we need an estimate of the temperature, density, or 
a combination. We can place the [Cnl/FIR ratio o f ~ .011 
derived in Sect. 3.2 on Fig. 11 in iBakes & Tielensl (1 19941) and 
find that the product Go^ff/n e « 20000. For a neutral re- 
gion, the electrons will come from the singly ionized carbon 
so n e ~ n H lA x 10~ 5 (see Sect. 3.3.3). Noting G 2 = G /100, 
T = 100r2, and n# = 1000«3, we obtain Gi y/T^/nj, « 1.5. 
In Sect. 3.1.2 we estimated Go ~ 50 b ut pointed out that this 
could be an underestimate by a factor 2. iKaufman et alj ( Il999t) 
estimate gas temperatures varying from 300 K to 100 K for 
G = 100 as densities climb from 100 to 10000 cm -3 , with 
tin = 1000 corresponding to T » 200 K. 

The general similarity between BCLMP 691 and 
BCLMP 302 holds also for the [Cnl/[Oi1 ratio, which is 
[Cn]/[0 1] ^ 1 in both regions (see Table 7 in iMookeriea et al.l 
1201 lj, and Nikola et al. in prep.). Looking at e.g. Fig. 
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4 of iKaufman et al.1 (Il999l) . we see that for a constant 
[Cii]/[Oi] ratio, decreases when T increases. The "solution" 
Go ~ 100, «h ~ 1000 seems to be a good fit and the value 
of T ~ 200K fits both with the IKaufman et all (1 19991) PDR 
model and with the [C n]/FIR ratio via the calculations in 
iBakes & Tielensl d!994l) . Looking back at Fig. H the [Cn] 
column density is about Afcn] ~ 1.3 x 10 17 ctrT 2 , in good 
agreement with Sect. 3. 3.2. It may be wor th noting that some 
"extra" H2 was found bv lBraine et alj d2010l) through correlating 
the Herschel/SPIRE cool dust measurements with CO and Hi 
observations. Their interpretation was indeed that this column 
density could likely be attributed to regions, cloud envelopes, 
where the CO was photo-dissociated but the gas remained 
predominantly molecular. The 1.3 ± xl0 17 cnT 2 C + column 
we observe comes from this layer of H 1 and H2 exterior to the 
region where most carbon is in CO. As the C + layer contains 
more H than is observed in Hi, particularly when considering 
only the Hi within the velocity range detected in [Cn], some 
molecular gas not detected via CO must be present. 

The correspondence along the cross between the [Cn] and 
CO intensities and velocities is quite good (see e.g. Fig. 0) in 
BCLMP 691 but this is not the case for all regions. We suspect 
that this may be a question of geometry and that in BCLMP 691 
we are probably seeing the H 11 region with the molecular cloud 
right behind it such that the surface of the cloud is in the same 
direction as the bulk of the cloud (where carbon is in CO). At 
the scale of this H 11 region, the [C n] emission is clearly com- 
ing from the PDR, plus perhaps the ionized skin. BCLMP 691 
appear s to be part of molecular cloud 299 in the iGratier et alJ 
(1201 lh catalog. The strong emission in all wavebands tracing 
star formation, including the FUV, suggests that indeed the H n 
region is on the near side of the GMC. While we use the term 
"cloud", molecular clouds have substructure so the exposed ma- 
terial may have a rather intricate perimeter. The CO peak is ac- 
tually just south-east of the Ha peak (the H 11 region) but the 
FUV, 8, and 24 //m emission all peak on the H 11 region, imply- 
ing that the warm part of the molecular cloud is behind the H n 
region. Since the [C 11] emission is coming from the near side of 
the molecular cloud, it is not surprising that it shares the same 
dynamics (velocities, Fig . [6j> . Perhaps due to the geometry, these 
observations of BCLMP 69 1 show clearly that the [C 11] emission 
is coming from the molecular cloud - Hn region interface, as 
in models of photo-dissoci ation regions dTielens & Hollenbachl 
1 19851: IKaufman et al.|[l999h . 
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Fig. 2. [C n] spectra superposed on a Spitzer 24 fim image at 6 " resolution of BCLMP 691. The dotted blue c ontours show the Ha 
emission at levels of 100, 300, 1000, 3000, and 10000 counts (lGreenawaitHl998tlHoopes & Walterbosll2000t) . Only positions with 
detected [C n] flux are shown. The velocity scale is from -300 to -230 km s _1 and the y axis is from -0.2 K to 1 K. The beamsize 
is shown as a circle around the southernmost spectrum. 
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Table 1. Positions observed in [Cn]. Columns 1 and 2 give the offset from the nominal pointing center of (J2000) RA: 01:34:16.40 
DEC: 30:51:54.6 in arcseconds (B1950 coordinates are RA: 01:31:26.98 DEC: 30:36:33.6). Columns 3, 4, and 5 give [Cn] fluxes, 
velocities and line widths for the detected positions. For the non-detections, column 3 gives the rms noise level in mK. CO integrated 
intensities, velocities, and line widths are provided in columns 6-8 and the H i column density is indicated in column 9. Velocities 
and line widths (and their uncertainties) are estimated via gaussian fits to the spectra; when no value is present, the routine did not 
converge or t he spe ctrum was not considered a det ection. The CO integra ted intensities and H i column densities are taken from 
iGratier et alJ (120101) . Finally, the 24 //m flux density (Tabatabaei et al. 2007) is given in the last column. 
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Fig. 7. Position-velocity plots along the two arms of the observed HIFI cross. Left: perpendicular to the major axis of the galaxy, 
Right: along the major axis of the galaxy. For each plot the color bitmap is the [C n] in main beam temperature (Kelvins), the white 
contours are the CO(2-l) Tmb temperature every 0.1 K starting from 0.1 K, the dark blue contours are the Hi temperature every 
10 K starting from 10 K. The horizontal bar at the bottom of each panel indicates the 12" beamsize of the observations. 



